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Abstract 


Based on the systematic investigation of the data available for A > 40, a 
Z*/3 dependence for the nuclear charge radii is shown to be superior to the 
generally accepted A!/? law. A delicate scattering of data around Rel Z 1/3 
is infered as owing to the isospin effect and a linear dependence of Re/Z 1/3 
on N/Z ( or (N — Z)/2 ) is found. This inference is well supported by the 
microscopic Relativistic Continuum Hartree-Bogoliubov (RCHB) calculation 
conducted for the proton magic Ca, Ni, Zr, Sn and Pb isotopes including the 
exotic nuclei close to the neutron drip line. With the linear isospin dependence 
provided by the data and RCHB theory, a new isospin dependent Z 1⁄3 formula 


for the nuclear charge radii is proposed. 
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Nuclear radius is one of the most fundamental bulk properties of an atomic nucleus 
[1,2]. Among all the size quantities describing nucleus, nuclear charge radii has been inves- 
tigated by various techniques and methods experimentally [3-9], including the muonic atom 
spectroscopy [3], isotope shift of optical and K X-ray spectroscopy [4-8] and high energy 
elastic electron scattering [9], etc.. Recently more and more nuclei far from the {-stability 
line become accessible experimentally thanks to the development of radioactive ion beam 
facilities [10,11]. The nuclear size connected with exotic phenomena such as skin and halo 
have become a hot topic. The understanding of its property has importance not only in 
nuclear physics, but also in other fields such as astrophysics and atomic physics, etc. With 
its accuracy, the study of the nuclear charge radii is very important to understand not only 
the proton distribution inside the nucleus but also the halo and skin. Particularly if one can 
get a simple and reliable formula for nuclear charge radii, it will be very useful to extract 
the de-coupling of proton and neutron in the exotic nuclei and provide information for the 
effective nucleon-nucleon interaction widely used in all the nuclear models. Here in this let- 
ter the available experimental charge radii data for A > 40 will be examined and its global 
behavior will be studied. Instead of the widely accepted Al’? law, a new Z'/3 formula with 
isospin effect will be proposed. 

Based on the consideration of the nuclear saturation property, nuclear charge radii Re 


are usually described by the A!/* law [1,2] 
Re = ra A, (1) 


where A is the mass number and R, = "a (r?) 1/2, with (r?)!/? the root-mean square (rms) 
charge radius. For very light nuclei, because of their small A and large fluctuation in charge 
distribution due to the shell effect with short period, it seems that the charge distribution 
radius as a bulk property has little meaning. A detail analysis of charge radius data for 
A > 40 shows that r4 is by no means a constant, but systematically decreases with A; i.e., 
ra © 1.31 fm for light nuclei (A ~ 40) and r4 = 1.20 fm for very heavy nuclei (see upper 


left panel in Fig. 1). This fact implies that some physics is missing in Eq. (1). 
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A definite evidence of the violation of A‘/? law is also found in the measurements of 
isotope shift in mean square charge radii [12,13]. In particular, 6(r?) 442,4 values (associated 
with an addition of two neutrons) are often found to be considerably smaller compared to 
what is expected from the A!/? law (Or?) Avo = sar?) a): A typical example is that the 
observed charge radii of the calcium isotopes “°~°°Ca remain almost the same (except a 
very little change induced by deformation or shell effect), though the mass number A has 
changed significantly. In contrast, there is also evidence that the observed 6(r?) 419.4 values 
(associated with the addition of two protons) are often greater than what is expected from 
the AY? law (e.g., d(r?) for “Ti—“Ca, »°Ti—Ca, etc.). 

Along the -stability line, the ratio Z/A gradually decreases with A, i.e., for light 
nuclei Z/A ~ 1/2, and for the heaviest -stable nucleus 738U, (Z/A)/> = 0.7285, thus 
(1/2)'/3/(Z/A)'/3 = 1.09, which is very close to the r4 ratio 1.30/1.20 shown in upper left 
panel of Fig. 1. A naive point of view is that the charge radius of a nucleus may be more di- 
rectly related to its charge number Z, rather than its mass number A. Therefore, compared 


to the Al/? law, a Z!/3 dependence for nuclear charge radii may be more reasonable 
R, = 7zZ"8. (2) 


as noted in Ref. [12]. An analysis of the very limited data of charge radii then available 
showed that rz remains almost a constant, i.e., rg = 1.65(2) fm for A > 40. The Z! 
dependence of nuclear charge radii was also used to modify the Coulomb energy term in 
the semi-empirical nuclear mass formula [14], and it was found that the agreement between 
the calculated and experimental results was improved. Moreover, the A~!/? law for the 
nuclear giant ( monopole, dipole and quadrupole ) resonance energy (x 1/R) also could be 
improved, if the A~!/? dependence is replaced by a Z~'/3 dependence [15]. 

In the past two decades, a vast amount of new experimental information on the elec- 
tromagnetic structure of nuclear ground states of many nuclei has become available [3-9], 
and accuracy has been improved. In particular the muon factories at Los Alamos (LAMPF) 


and at Villigen (PSI, formerly SIN) started their operation at 1974. Almost all stable nuclei 
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have been measured by the muonic X-ray transition technique and the corresponding charge 
radii have been rather accurately deduced (the experimental relative error is about 107°). 
Moreover, modern techniques for optical isotope shift measurements have made it possible 
to reach even short-lived (down to 1 s) unstable isotopes [3]. Therefore, it is worthwhile to 
reexamine the fundamental property of nuclei and to investigate whether the vast amount of 
improved experimental results follow the Z!/3 dependence. The values of measured (r?)1/? 
for 536 nuclei with A > 40 compiled in Ref. [3-9] are analyzed in Fig. 1 by using the 
A'/3 and Z‘/3 dependence, respectively. The dependence of charge radii on the quadrupole 


deformation @ has been taken into account for the rare-earth deformed nuclei, i.e. [1] 
ra=Tadit > 62) rg =fząa(l + > 62) (3) 
A=TAd a h Z= TZzd m 


and for spherical nuclei (8 = 0): r4 = raga, Tz = Tza, and the values of 8 are taken from 
Ref. [9,16]. 

In the upper left and right panels of Fig.1, the charge radii for the most stable 159 
nuclei with A > 40 along the 6-stability line have been analyzed by using the AV’? and Z1/3 
dependence. In the middle left and right panels, the same has been done for the measured 
(r?\1/2 for 536 nuclei with A > 40. Two significant features can be observed: (A) On the one 
hand, the agreement between the data and the calculated results using the Z'/? dependence 
is much better than that using the A!/? law, i.e., while there exists a global regular decrease 
of raa with A, rza nearly remains constant (rzą = 1.631(11) fm). The relative rms deviations 
o for the Z'/3 dependence (o = 7.57 x 107° for stable nuclei and 1.00 x 107? for 536 nuclei 
) are much less than that for the A'/? law (o = 1.90 x 107? for stable nuclei and 1.63 x 107? 
for 536 nuclei ). (B) On the other hand, though the rms deviation for the Z'/* dependence 
is significantly reduced, an isospin induced scattering of the data in the middle panels in 
Fig. 1 can be also observed compared with that in the top panels. In fact, rzą generally 
increases with N for most isotopic chains, e.g. for 4) °°Zr, 23-'°°Mo, #8-"*Ru, 30?-1°Pd, 
106-116 Qyq 112—124g), 122-130, 124-136%e 2-18 144-1545m, 154-160Gq etc, (except for 


only a few lighter isotopic chains, e.g., 8 °°Kr, §3-°8Sr, and a small anomalous decrease of 


rza with N due to the shell closure at N = 50 is observed). Therefore, it seems necessary 
to investigate an isospin dependent correction for the scattering of rza. In Ref. [13,17], the 
isospin effect has been considered based on the A!/? law. However, considering the fact that 
the Z'/3 dependence can describe the nuclear charge radii much better than the A! law, 
we take the Z'/3 dependence as a more reasonable starting point for describing the isospin 
dependence of nuclear charge radii. 

To confirm that isospin dependent Z! formula to be developed for nuclear charge radii 
mentioned above is also valid for nuclei far from the {-stability line, the charge radii in 
exotic nuclei is needed. However, as no such data is available, what we can do is to require 
that our new isospin dependent Z'/? formula should assort with a reliable and microscopic 
nuclear model. 

The fully self-consistent and microscopic relativistic continuum Hartree-Bogoliubov 
(RCHB) theory, which is an extension of the relativistic mean field (RMF) [18-20] and 
the Bogoliubov transformation in the coordinate representation [21], is a good candidate 
for the present purpose. The RCHB theory can describe satisfactorily the ground state 
properties for nuclei both near and far from the -stability line. A remarkable success of 
the RCHB theory is the self-consistent reproduction of the halo in !Li [22] and prediction 
of the exotic phenomenon - giant halo [23]. In combination with the Glauber model, the 
RCHB theory successfully reproduces the interaction cross section in Na isotopes [24] and 
the charge changing cross section of C, N, O, F isotopes (ranging from the (-stability line to 
the neutron drip line) on the target of !*C at 930 MeV/u [25,26]. These successes encourage 
us to apply the RCHB theory for the description of charge radii of nuclei both close to 
and far from the 6-stability line and check its validity for the data available and provide 
information for nuclei far away from the stability line. 

The detailed formalism and numerical techniques of the RCHB theory can be found in 
Ref. [21] and the references therein. In the present calculations, we follow the procedures 
in Ref. [21,23,24] and solve the RCHB equations in a box with the size R = 20 fm and a 


step size of 0.1 fm. The parameter set NL-SH [27] is used, which aims at describing both 


K 


the stable and exotic nuclei. The density dependent 6-force in the pairing channel with 
po = 0.152 fm~ is used and its strength Vo is fixed by the Gogny force as in Ref. [21]. The 
contribution from continua is restricted within a cut-off energy Eeuw ~ 120MeV. 

As typical examples, we studied the even-even Ca, Ni, Zr, Sn and Pb isotopes ranging 
from the -stability line to neutron drip line. The two neutron separation energies Sən 
is one of the essential quantities to test a nuclear model. In Fig. 2, the calculated Sən 
(open symbols) of the even-even Ca, Ni, Zr, Sn and Pb isotopes by the RCHB theory are 
compared with the data available (solid symbols) [28], where a satisfactory agreement is seen. 
Particularly the deviation between the calculated binding energies with the data available 
is within 1%. In the present calculation, the neutron drip-line nuclei are predicted at ‘Ca, 
INi, °Zr and Sn, respectively. In the Sən versus N curve for each isotopic chain, there 
are some kinks due to the neutron shell or subshell closure. For example, the closed shells 
at N = 20, 28 and subshell at N = 40 correspond to kinks in the S», versus N curve for 
Ca isotopes at “Ca, *8Ca and Ca, respectively. While the kink at N = 20 for Ca may 
be also due to the Wigner term for N = Z = 20. However, there are no kinks at "Ca and 
176Sn, which indicate the disappearance of magic number N = 50 and 126 for these nuclei 
in RCHB. 

The rms charge radii (r?)!/? obtained from the RCHB theory (open symbols) and the 
data available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb isotopes are given in 
Fig. 3. As it could be seen, the RCHB calculations reproduce the data very well ( within 
1.5%). For a given isotopic chain, an approximate linear N dependence of the calculated 
rms charge radii (r?)!/? is clearly seen in Fig. 3, which shows that the variation of (r?)!/? 
for a given isotopic chain deviates from both the the simple Z!/* dependence and the simple 
A! law (denoted by dashed lines in Fig. 3). Therefore, a strong isospin dependence of 
nuclear charge radii is necessary for nuclei with extreme N/Z ratio. 

In Fig. 4, the experimental and RCHB predicted rzą = R./Z'/? for the proton magic 
isotopes are presented as a function of n = N/Z. It is clearly seen that the coefficient rza 


increases linearly with 7 (except some deviations due to deformation or shell effect) and the 


A 


slopes are nearly the same for these isotopic chains. The linear y ( or isospin Tz = (N—Z)/2 
) dependence of rzq for an isotopic chain may be understood as the effect of the first order 
perturbation correction of nuclear wave function due to an isospin Tz dependent interaction 
[29]. Based on the analysis of data in the middle and upper panels of Fig. 1 and RCHB 
prediction in Fig. 3 and 4, we propose the following isospin dependent Z!/? formula for 


nuclear charge radii: 
Re=a ZP +b- n), n= N/Z, (4) 


where 7* is 7 = N/Z for the nuclei along the (6-stability line which can be directly extracted 
from the nuclear mass formula [1], a = r34(1 + 267), rz4 = 1.631(11) fm obtained in upper 
right panel of Fig.1, and b = 0.062(9) obtained from the least square fitting. 

The analysis of the available data using Eq. (4) with r},, and b thus obtained is displayed 
in the lower right panel of Fig. 1. It is found the data are reproduced better by Eq. (4) than 
by Eq. (2) (the rms deviation is reduced by about 40%). The same has been done for Al? 
dependence but with less success ( see the lower left panel of Fig. 1 ). In Ref. [13,17], a similar 
equation for A!/? dependence has been used to describe the isospin dependence of charge 
radii and a better agreement has been achieved by fitting both rą and b simultaneously. 
But then a simple explanation for r4 (saturation property) and b(isospin effect) in Eq. (4) 
is missing. It is expected that the modified Z! formula (Eq. (4)) will become more useful 
with more and more data obtained for the nuclei far from the (-stability line. 

In summary, we have systematically investigated the nuclear charge radii with A > 
40. It is clearly seen that the Z'/3 dependence is superior to the A? law. A delicate 
scattering of data around R,/Z'/3 = 1.631 is infered as owing to the isospin effect and a 
linear dependence of R,/Z'/3 on N/Z ( or (N — Z)/2 ) is found. This inference are well 
supported by the microscopic RCHB calculation conducted for the proton magic Ca, Ni, 
Zr, Sn and Pb isotopes including the exotic nuclei close to the neutron drip line, which 
reproduce Sən and nuclear charge radii data available well. With the linear dependence 


of the coefficient rzą on N/Z (or (N — Z)/2 ) read from data and RCHB theory, a new 
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isospin dependent Z!/° formula for the nuclear charge radii is proposed, which improves the 
description of the data available for nuclei near the (-stability line and could be very useful 
for new data obtained for nuclei far from the 6-stability line. 
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FIG. 1. The nuclear charge radius data for r4 in A! and rz in Z'/3 law with and without 


isospin dependence, for the details see the text. 
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FIG. 2. Two-neutron separation energies Son of even Ca, Ni, Zr, Sn, Pb isotopes as a function 


of N, including the data (solid symbols) from Ref.[28] and the RCHB calculation with a 6-force 


(open symbols). 
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FIG. 3. The rms charge radii versus the neutron number N for even-even Ca, Ni, Zr, Sn, Pb 
isotopes. The RCHB calculation with 6-force is represented by open symbols, while the corre- 


sponding data is denoted by solid symbols. The dashed lines represent the predictions by the A!/? 


law with rag = 1.228 fm. 
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FIG. 4. The experimental (solid symbols) and RCHB predicted ( dashed lines ) coefficient 
rzą = R¿/Z™ for the nuclear charge radii as a function of isospin quantity 7 = N/Z in even-even 


Ca, Ni, Zr, Sn and Pb isotopes. An asymptotic behavior is drawn as a solid line. 


